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L  Introduction 


In  systems  where  independent  users  transmit  through  a  single  common  channel,  the  deploy¬ 
ment  of  random  access  transmission  algorithms  is  frequently  desirable,  for  the  following  rea¬ 
sons:  (1)  They  are  implemented  independently  by  each  user,  without  a  prior  coordination 
among  the  users.  (2)  They  are  insensitive  to  changing  user  population.  (3)  They  induce  low 
delays  wfcen  the  user  traffic  is  bursty. 

In  this  paper,  we  present  and  analyze  a  full  sensing  window  random  access  algorithm.  The 
algorithm  was  first  proposed  for  systems  with  strict  delay  limitations,  [3],  and  it  requires  that 
each  user  know  the  overall  feedback  history,  (full  sensing).  As  compared  to  other  such  existing 
algorithms,  the  present  algorithm  has  the  following  interesting  and  beneficial  properties:  (1)  It 
can  be  easily  modified  to  operate  in  limited  sensing  environments,  where  each  user  follows  the 
feedback  history  from  the  time  he  generates  a  message  to  the  time  when  this  message  is  success¬ 
fully  transmitted.  (2)  In  the  presence  of  the  limit  Poisson  user  model,  the  algorithm  attains  the 
same  throughput  as  that  attained  by  the  Capetanakis’  dynamic  algorithm,  [1],  while  it  induces 
significantly  lower  delays  for  arrival  rates  above  0.30,  and  superior  resistance  to  feedback  errors. 
(3)  The  simple  operations  of  the  algorithm  allow  analysis  and  evaluation  when  strict  delay  limi¬ 
tations  exist,  [3].  Its  simplicity,  in  conjunction  with  its  regenerative  properties,  provide  the 
means  for  the  analysis  and  evaluation  of  the  output  traffic  interdeparture  distribution  induced  by 
the  algorithm.  The  analysis  of  the  latter  distribution  is  important  when  several  systems  which 
use  some  Random  Access  Algorithm,  (RAA),  for  internal  transmissions  interact,  and  it  is  not 
quite  feasible  when  either  the  Capetanakis,  [1],  or  the  Gallager,  [2],  algorithms  are  deployed.  (4) 
As  compared  to  Gallager’s  algorithm,  [2],  the  present  algorithm  operates  in  environments  where 
the  Poisson  model  is  not  valid,  (e.g.,  when  more  than  one  packets  can  be  generated  within  a 
given  time  instant),  and  can  be  then  analyzed. 

The  organization  of  the  paper  is  as  follows:  In  section  II,  the  system  model  is  presented, 
and  the  algorithm  is  described.  In  section  III  the  throughput  and  delay  analyses  arc  included,  in 
the  presence  of  the  limit  Poisson  user  model,  and  in  the  absence  of  feedback  errors.  In  section 
IV,  the  performance  of  the  algorithm  in  the  presence  of  feedback  errors  and  its  operations  in  lim¬ 
ited  sensing  environments  are  discussed.  In  section  V,  the  output  traffic  interdeparture  distribu¬ 
tion  induced  by  the  algorithm  is  analyzed  and  evaluated.  In  section  VI,  some  conclusions  are 
drawn. 

n.  The  System  Model  and  the  Algorithm 

We  assume  packet  transmitting  users,  slotted  channel,  binary  collision  versus  noncollision, 
(CNC),  feedback  per  slot,  no  propagation  delays,  and  absence  of  feedback  errors.  We  also 
assume  nonexistence  of  error  correction  coding;  thus,  collided  packets  are  fully  destroyed  and 
retransmission  is  then  necessary.  Time  is  measured  in  slot  units,  slot  t  occupies  the  time  interval 
[t,  t+1),  and  x,  denotes  the  feedback  that  corresponds  to  slot  t;  xt  =  C  and  xt  =  NC  represent  then 
collision  and  noncollision  slot  t,  respectively.  For  this  system,  let  the  following  full  sensing  syn¬ 
chronous  random  access  algorithm  be  deployed. 

The  algorithm  utilizes  a  window  of  length  A.  Let  t  be  a  time  instant  such  that,  for  some 
t,  <  t,  all  the  packet  arrivals  in  (0,  tj  have  been  successfully  transmitted  by  the  algorithm  and 
there  is  no  information  regarding  the  arrival  interval  (tj,  t],  and  such  that  t  corresponds  to  the 


beginning  of  some  slot.  The  instant  t  is  then  called  Collision  Resolution  Point,  (CRP),  the 
arrival  interval  (0,  tj]  is  called  "resolved  interval",  and  the  interval  (tj  t]  is  called  "the  lag  at  t". 

In  slot  t,  the  packet  arrivals  in  (tit  t2  ^  min  (t!  +  A.t)]  attempt  transmission,  and  the  arrival  inter¬ 
val  (tj,  tj  is  then  called  the  "examined  interval".  If  (tj,  t2]  contains  at  most  one  packet,  then  it  is 
resolved  at  t.  If  (tj,  t2]  contains  at  least  two  packets,  instead,  then  xt  =  C,  a  collision  occurs  at  t, 
and  its  resolution  starts  with  slot  t+1.  Until  the  collision  at  t  is  resolved,  no  arrivals  in  (t2,  are 
allowed  transmission.  The  time  period  required  for  the  resolution  of  the  latter  collision  is  called 
the  Collision  Resolution  Interval,  (CRT).  During  some  CRI,  each  user  acts  independently  via  the 
utilization  of  a  counter  whose  value  at  time  t  is  denoted  rt.  The  counter  values  can  be  either  1  or 
2,  and  they  are  updated  and  utilized  according  to  the  rules  below. 

1.  The  user  transmits  in  slot  t,  if  and  only  if  rt=  1.  A  packet  is  successfully  transmitted  in  t,  if 
and  only  if  r,  =  1  and  xt  =  NC. 

2.  The  counter  values  transition  in  time  as  follows: 

(a)  If  xt_i  =  NC  and  rt_t  =  2,  then  rt  =  1 

(b)  If  xt_[  =  C  and  r,_1  =  2,  then  rt=2 

(c)  If  x,_j  =  C  and  r,_j  =  1,  then 


«■«= 


1. with  probability  0.5 

2, with  probability  0.5 


Remarks  We  note  that  the  algorithmic  operations  can  be  depicted  by  a  two-cell  stack,  where  at 
each  time  instant  t,  ceil  1  contains  the  transmitting  users,  (those  with  rt  =  1),  and  cell  2  contains 
the  withholding  users,  (those  with  rt  =  2).  The  algorithm  lumps,  thus,  the  unsuccessful  users 
together.  In  contrast,  Capetanakis'  algorithm  distributes  the  unsuccessful  users  across  the  cells 
of  an  infinite-cell  stack.  As  with  Capetanakis’  dynamic  algorithm,  the  window  size  A  is  here 
subject  to  optimization  for  throughput  maximization. 


HI.  Throughput  and  Delay  Analysis 

In  this  section,  we  present  the  throughput  and  delay  analyses  of  the  algorithm,  in  the 
absence  of  feedback  errors,  in  the  full  sensing  environment,  and  in  the  presence  of  the  limit  Pois¬ 
son  user  model.  As  proven  in  [51,  the  latter  user  model  represents  a  lower  bound.  That  is,  when 
the  user  population  in  finite,  the  users  are  independent  and  identical,  and  the  packet  generation 
process  per  user  is  memoryless,  then  the  throughput  and  delay  characteristics  of  the  algorithm 
are  superior  to  those  induced  when  the  user  environment  is  limit  Poisson. 

Let  the  system  start  operating  at  time  zero,  and  let  us  consider  the  sequence  in  time  of  lags 
that  are  induced  by  the  algorithm.  Let  Q  denote  the  length  of  the  i-th  lag,  where  i>l.  Then,  the 
first  lag  corresponds  to  the  empty  slot  zero;  thus,  Ct  =  1.  In  addition,  the  sequence  C,;  i>l  is  a 
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•O 


’i 
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Markov  chain  whose  state  space  is  at  most  countable.  Let  Dn  denote  the  delay  experienced  by 
the  n-th  successfully  transmitted  packet  arrival,  as  induced  by  the  algorithm;  that  is,  the  time 
between  the  arrival  of  the  packet  and  its  successful  transmission.  Let  the  sequence  Tj,  i>l  be 
defined  as  follows:  Each  Tj  corresponds  to  the  beginning  of  some  slot,  and  T,  =  1.  Also,  each  Tj 
corresponds  to  the  ending  point  of  a  length-one  lag.  is  then  the  ending  point  of  the  first  after 
Tj  unity  length  lag.  Let  Rj,  tel  denote  the  number  of  successfully  transmitted  packets  in  the  time 
interval  (0,  TJ.  Then,  Qj  =  Ri+1  -  Rj,  tel  denotes  the  number  of  successfully  transmitted  packets 
in  the  interval  (Tj,  Tj+1].  The  sequence  Q,,  tel  is  a  sequence  of  i.i.d.  random  variables;  thus 
Rj,  tel  is  a  renewal  process.  In  addition,  the  delay  process  Dn,  rtel  induced  by  the  algorithm  is 
regenerative  with  respect  to  the  process  Rit  tel  and  the  distribution  of  Qj  is  nonperiodic,  since 
P  (Qi  =  1)  >  0. 

Let  us  define, 

Z=E{Q1},  W  =  E{ ]£  Dj }  (1) 

i=l 

From  the  regenerative  arguments  in  [4],  it  follows  that  the  expected  per  successfully 
transmitted  packet  steady-state  delay,  D ,  is  given  by  the  following  expression: 

D  =  W271  (2) 

The  effective  computation  of  D  relies  on  the  successful  derivation  of  upper  and  lower  bounds  on 
the  quantities  W  and  Z.  Those  bounds  are  found  via  the  utilization  of  the  methodology  in  [4],  in 
conjunction  with  the  quantities  defined  in  [3];  the  details  are  thus  omitted.  If  E  { / 1  A,d }  denotes 
the  expected  length  of  a  CRI,  given  that  it  starts  with  an  examined  interval  of  length  A  and  with 
a  lag  d,  then  bounds  on  W  and  Z  can  be  found  only  if: 

A>E{/IA,d)  (3) 

The  inequality  in  (3)  determines  the  stability  region  of  the  algorithm.  Let  us  define, 

Ln.k-n:  The  expected  number  of  slots  needed  by  the  algorithm  for  the  successful 

transmission  of  k  packets,  given  that  n  of  the  k  packets  have  counter  values  equal 
to  1  and  that  the  remaining  k-n  packets  have  counter  values  equal  to  2. 


Then,  as  found  in  [3],  we  have  the  following  expressions,  where  X  denotes  the  intensity  of  the 
Poisson  traffic  process: 


0  <  Lv  o  £  —  k2+-7*k-2  ,  k  >  1 

(4) 

Lo.o  =  =  1  »  L0ii  =  1  +  Lj  o  ,  tel 

(5) 

E{/IA,d}=SLki0e-^-^|^ 

k=0  K • 

(6) 

The  expressions  in  (4),  (5),  and  (6),  in  conjunction  with  the  methodology  in  [4],  are  used  in  the 


computation  of  the  algorithmic  throughput  and  delays  when  the  limit  Poisson  user  model  is 
present.  In  Table  1,  we  include  the  computed  upper  and  lower  bounds,  Du  and  D;  respectively, 
on  the  expected  per  packet  delay  D,  for  various  Poisson  intensities  X,  and  for  both  the  present 
and  the  Capetanakis  dynamic  algorithms. 


Proposed  a 

lgorithm 

Capetanakis  dynamic 
algorithm 

X 

D' 

Du 

Dl 

Du 

0.02 

1.562 

1.563 

1.563 

1.564 

0.06 

1.708 

1.716 

1.713 

1.719 

0.10 

1.888 

1.917 

1.903 

1.921 

0.16 

2.257 

2.363 

2.308 

2.362 

0.20 

2.607 

2.812 

2.712 

2.809 

0.24 

3.103 

3.467 

3.308 

3.476 

0.30 

4.412 

5.197 

4.976 

5.365 

0.32 

5.162 

6.170 

5.973 

6.501 

0.36 

7.941 

9.665 

9.798 

10.883 

0.38 

11.008 

13.398 

14.121 

15.855 

0.40 

18.262 

22.024 

24.427 

27.736 

0.42 

57.354 

67.665 

78.530 

90.212 

Table  1 


Upper  and  Lower  Bounds  on  Steady-State  Expected  Delays 


Regarding  the  throughput  X*  and  the  optimal  window  size  A*,  the  following  results  were 
found. 


_ Proposed  Algorithm: _ X  =  0.4295  A  =  2.33 

Capetanakis*  Dynamic  Algorithm:  X *  =  0.4295  A*  =  2.673 

Table  2 

Throughputs  and  Optimal  Window  Sizes 

From  Table  2,  we  observe  that,  in  the  presence  of  the  limit  Poisson  user  model,  the  algo¬ 
rithm  in  this  paper  attains  the  same  throughput  with  the  Capetanakis  dynamic  algorithm,  but  util¬ 
izes  a  smaller  window  size.  From  Table  1,  we  observe  that  the  two  algorithms  induce  practically 
identical  delays  for  Poisson  rates  in  (0,  0.30),  while  for  Poisson  rates  in  (0.30,  0.42],  the  present 
algorithm  is  significantly  superior. 

Remarks  It  may  seem  surprising  that  the  algorithm  in  this  paper  attains  the  same  throughput  with 
the  Capetanakis  algorithm,  and  that  it  outperforms  the  latter  in  terms  of  delay  performance. 
Indeed,  the  expected  lengths  Lt0  in  (4)  are  bounded  by  quadratic  expressions,  while  the  same 
lengths,  Lk,  in  the  Capetanakis  algorithm  are  bounded  by  linear  functions  of  k.  However, 
L^o  =  4.5  while  L2  =  5.  At  the  same  time,  since  A*  =  2.33  for  the  present  algorithm,  the  proba¬ 
bility  of  a  higher  than  two  multiplicity  collision  is  very  small.  The  multiplicity-two  events  thus 
prevail,  and  the  algorithm  in  this  paper  becomes  superior  to  the  Capetanakis’  algorithm.  We 
note  that,  as  found  in  [3],  the  algorithm  performs  very  well  in  environments  where  strict  delay 
limitations  exist.  Then,  it  allows  significant  improvement  in  delay  performance,  at  the  expense 
of  minimal  loss  in  traffic.  In  addition,  the  analysis  of  the  algorithmic  performance  when  strict 
delay  limitations  exist  is  relatively  simple,  while  the  same  analysis  for  the  algorithms  in  [1]  and 
[2]  is  then  exceedingly  complex. 


IV.  Performance  in  the  Presence  of  Feedback  Errors  and  Operations  in  the  Limited 
Sensing  Environment 

Due  to  its  simple  rules,  the  present  algorithm  is  much  less  sensitive  to  feedback  errors  than 
the  Capetanakis’  algorithm.  To  see  that,  let  8  >  0  be  the  probability  that  an  empty  slot  is  errone¬ 
ously  interpreted  as  a  collision  slot,  and  let  this  be  the  only  form  of  feedback  error.  Then,  as 
shown  in  [7],  the  throughput  of  the  Capetanakis’  dynamic  algorithm  reduces  to  zero,  if  5  exceeds 
0.5.  In  contrast,  the  throughput  Xg  of  the  present  algorithm  remains  positive,  for  every  8  in  [0,  1 
).  Indeed,  given  8,  let  Bk,  Ek ,  and  Sk  respectively  denote  the  expected  number  of  collision, 
empty,  and  success  slots  during  the  resolution  of  a  k-multiplicity  collision,  and  let  Lk  be  the 
expected  number  of  slots  needed  for  the  resolution  of  this  collision.  Then,  the  simple  operations 
of  the  algorithm  easily  induce  the  following  expressions: 

Bk  =  Bk  ,  Sk8  =  Sk°  ,  Ek  =(1-S)-1  Ek  ;  0  <  8  <  1  ,  k  >  2 


L®  =  B®  +  Sk°  +  Ek  (1-8)-1  =L^  +  Sd-Sr1  Ek° £ (1-8)"1  Lk°;0^5<l,k>2 

Lo  =  (1-5)-2  ;  0  £  5<  1 

X8*  =  xTo(xf5_1  (x)) 


where 


uwi  £48e'xrr  •  4s =L 


Thus, 


.  X5*  Xo*  (1-8)2  =  (1-5)2  (0.4295)  >0  ;  for  all  8  in  [0,1) 


Using  methods  as  in  [3],  we  developed  bounds  on  the  expected  lengths  L®,  and  then  bounds 
on  the  throughput  Xg,  for  8e(0, 1).  We  applied  similar  methods  to  compute  bounds  on  the  same 
throughput  for  the  Capetanakis  dynamic  algorithm,  (in  [7]  the  Capetanakis  nondynamic  algo¬ 
rithm  is  considered).  We  exhibit  our  results  in  Figure  1.  We  obsrve  the  uniform  superiority  of 
the  proposed  algorithm.  We  notice  in  particular,  that  for  8  =  0.5,  the  proposed  algorithm  attains 
throughput  0.325,  while  the  throughput  of  Capetanakis’  dynamic  algorithm  is  then  zero. 

We  next  considered  the  case  when  a  single  transmission  may  be  interpreted  as  a  collision, 
with  probability  e,  and  when  this  is  the  only  type  of  feedback  error.  Using  the  same  methodol¬ 
ogy  as  above,  we  then  computed  bounds  on  the  induced  throughput,  X*,  for  both  the  proposed 
and  the  Capetanakis  dynamic  algorithms.  Our  results  are  included  in  Table  3.  We  observe  that 
the  two  algorithms  have  practically  identical  performance  in  this  case. 


*  r  “ 


Proposed  algorithm 

Cape  tan  akis  dynamic  algorithm 

e 

_ £ _ 

K 

0.02 

| 

0.422 

0.04 

■  '  • 

0.414 

0.06 

■  ■„  VB 

0.407 

0.08 

0.395 

0.399 

0.10 

0.387 

0.392 

0.20 

0.346 

0.353 

0.30 

0.305 

0.313 

0.40 

0.265 

0.273 

0.50 

0.225 

0.231 

0.60 

0.184 

0.187 

0.70 

0.142 

0.142 

0.80 

0.098 

0.097 

0.90 

0.052 

0.049 

1.00 

0.000 

0.000 

I  Table  3 


Throughput  when  with  probability  e,  a  single  transmission  is  interpreted  as  a  collision. 

Operations  in  the  Limited  Sensing  Environment 

In  the  limited  sensing  environment,  it  is  required  that  each  user  monitor  the  channel  feed¬ 
back  only  from  the  time  he  generates  a  packet,  to  the  time  this  packet  is  successfully  transmitted. 
Therefore,  the  users’  knowledge  of  the  channel  feedback  history  is  then  asynchronous.  The 
objective  in  this  case  is  to  prevent  new  arrivals  from  interferring  with  some  collision  resolution 
in  progress.  This  is  possible,  if  each  user  can  decide  whether  a  collision  resolution  is  in  progress 
or  not,  within  a  finite  number  of  slots  from  the  time  he  generates  a  new  packet.  We  observe  that 
a  user  who  has  a  new  packet  and  observes  a  C  slot  decides  to  wait,  since  he  can  then  deduce  that 
there  is  some  collision  resolution  in  progress.  Also,  since  a  CRI  ends  with  two  consecutive  NC 
slots,  all  the  users  who  observe  such  an  event,  decide  that  there  is  no  collision  resolution  in  pro¬ 
gress.  In  view  of  the  above  observations,  we  conclude  that  in  the  limited  sensing  environment, 
the  algorithm  can  be  modified  to  operate  as  follows: 


The  window  size  is  the  same  as  in  the  full  sensing  case.  The  window  slides  from  present  to  past, 
however.  In  particular,  the  edge  of  the  window  is  maintained  one  slot  before  the  present  time, 
and  the  window  slides  through  the  unexamined  interval  from  present  to  past,  (see  Fig.  2). 
Within  each  window,  the  operations  of  the  algorithm  are  the  same  as  in  the  full  sensing  environ¬ 
ment. 


In  the  limited  sensing  environment,  and  for  very  light  input  traffic,  the  algorithm  induces 
expected  per  packet  delay  equal  to  2.5.  As  the  rate  of  the  input  traffic  increases,  the  expected 


delays  approach  those  induced  in  the  full  sensing  environment  The  throughput  ot  the  aa'*ir:"  - 
remains  identical  to  that  in  the  full  sensing  environment  In  Figure  L  we  plot  the  r«pr  : 
delays  that  the  algorithm  induces,  in  both  the  full  sensing  and  the  limited  sensing  ms  imnrn'-- 
In  the  latter  environment,  the  expected  delays  were  computed  via  methodologies  a'  :r.  >se  • 

(9). 


Remarks  We  point  out  that  the  modification  of  Capetanakis  dynamic  aig  -r.:-  •:  • 
in  the  limited  environment,  is  still  an  open  and  complex  problem  If  i-,::.. 
modification  is  simple  when  the  proposed  algorithm  is  adopted  lr.  sssicm  ±'x:- 
user  model  is  valid,  the  Pan -and  Try  algonthm  with  b.nars  feedback  r  »  ar  •» 
operate  in  the  limited  sensing  environment.  |9|.  |1(>)  The  throughput  >t  :r,e 
then  0.45.  But  when  the  Poisson  user  model  ’s  not  salid  it  leads  to  dead,  *.*, 
proposed  algonthm  performs  well  in  nonPotsson  user  systems 


V.  The  Output  Traffic  Interdeparture  Distribution 

In  this  section,  we  concentrate  on  the  computation  of  the  outpu.  traf:.*  .r.rer.:-; 
bution.  In  particular,  we  find  analytically  the  steady-state  distribution  >r  trie  .:.r  - 
two  consecutive,  successful  transmissions,  when  the  algonthm  in  this  paper  a 

point  out  that  the  algonthm  generates  an  output  traffic  process  with  memory  •  >ar  , 
correspond  thus  to  the  first  order  distribution  from  this  process  This  first  order  d»str. 
conjunction  with  a  memoryless  assumption,  can  be  used  as  an  approximation  of  rhe  .u\.a. 
traffic  process,  when  studies  of  systems  which  deploy  the  algonthm  and  interact  w  itr,  ea.* 
are  undertaken.  Such  interactions  may  correspond,  for  example,  to  ill  Servicing  the  >  .. 

traffic  from  several  systems  that  deploy  the  algonthm,  by  a  single  server  queue  2>  Trans 
ting  the  output  traffic  from  a  system  that  deploys  the  algonthm,  through  the  transmission  «.  r.ar.r.e. 
of  another  random  access  system,  (multi-hop  problem).  The  methodology  wc  use  to  compute  tree 
first  order  distribution  from  the  output  traffic  process,  extends  easily  to  higher  order  distribution^ 
from  the  process  as  well.  The  computations  become  then  exceedingly  complex,  however 

Our  methodology  utilizes  the  regenerative  character  of  the  output  traffic  process  that  the 
algorithm  generates,  and  its  steps  are  as  those  in  [4],  The  initial  challenge  here  lies  in  the  deter 
mination  of  regenerative  points,  which  are  pertinent  to  the  output  traffic  process.  We  define  the 
sequence  {P,}^  of  such  points  as  follows:  Each  Pj  is  a  collision  resolution  point,  (CRP),  w  hich 
follows  a  slot  containing  a  successful  transmission  and  at  which  the  lag  equals  one.  P,  is  the  first 
after  zero  such  CRP,  and  for  every  i>l,  Pi+1  is  the  first  after  P,  such  CRP.  Let  S,  i>l  denote  the 
number  of  successful  transmissions  in  (0,  PJ,  and  let  d,,  denote  the  distance  between  the  (n- 1  )-th 
and  the  n-th  successful  transmission.  Then,  S,,  i>I  is  a  renewal  process,  and,  as  it  can  be  easily 
seen,  the  process  d,,  n>l  is  regenerative  with  respect  to  it.  Let  us  define,  Cj  =  Si+1  -- S,,  i>l. 
Then  C;  denotes  the  number  of  successful  transmissions  in  the  internal  ( P,,  P1+1],  where  this 
interval  will  be  called  the  i-th  cycle.  Let  us  define, 

fl  i  if  dn  =  s 
In(s)  -  j  o  ,  otherwise 


(10) 


H  =  E(Pi+1-Pi} 


ifinwvi 


From  the  regenerative  theorem  [4],  we  then  conclude  that  if  C  =  E{Cj  }< then, 


lim  ,  n  lim  N  ,  Q 

N_*.  N-1  Z  E(  £  I„(s))  =cr*  E(  £  Us)) 

n=l  n=l  n=l 


(12) 


;  where  if  the  intensity  of  the  input  Poisson  traffic  is  X,  then, 

C  =  XH 


(13) 


In  addition,  since  P(Cj  =  1)  >  0,  the  distribution  of  is  aperiodic  and  there  exists  a  random 
variable  d„,  such  that  the  sequence  d,,,  n=l,2,...  converges  in  distribution  to  d^.  Then,  dM 
represents  the  steady  state  interdeparture  distance  induced  by  the  algorithm,  and  its  distribution 
satisfies  the  equality. 


* 


P(d0.  =  s)  =  C"1E{£  ^(s)) 


(14) 


n=l 


The  finiteness  and  the  computation  of  the  quantities  C  and  E{  £  ^(s)}  in  (14)  are  related  to  the 

n=l 

existence  and  computation  of  appropriate  solutions  to  infinite -dimensionality  linear  systems. 
Those  systems  and  their  solutions  are  included  in  the  Appendix.  In  Table  4,  we  include  the  com¬ 
puted  upper  and  lower  bounds,  respectively  denoted  P“  and  P on  the  probability  P(d00  =  s),  for 
input  traffic  Poisson  intensities  A.  =  0.1  and  X  =  0.4.  In  Figure  4,  we  plot  the  lower  bounds 
against  s,  for  various  input  traffic  Poisson  intensities  X. 
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X  =0. 1  and  X  =0.4 


From  the  results  in  Table  4  and  Figure  4,  we  draw  the  following  conclusions:  (1)  For  low 
rates  X  of  the  Poisson  input  traffic,  (X  S0.1),  the  interdeparture  distribution  is  close  to  the  Ber¬ 
noulli  distribution  whose  parameter  is  p  =  Xe-\  In  particular,  denoting  P,  ^  P(d„  =  s),  we  then 
have  Ps~p  (l-p)s-1,for  s  £  2.  The  probability  P1(  however,  is  then  significantly  larger  than  the 
Bernoulli  parameter  p.  The  intuitive  explanation  of  the  latter  phenomenon  goes  as  follows:  For 
small  rates  X,  single  arrivals  in  two  consecutive  slots  occur  with  probability  p2  =  (Xe~*)2  =  X2, 
while  the  probability  of  a  collision  slot  is  then  approximately  equal  to  2_1  X2  e_x  =  2“'  X2. 
Thus,  for  small  rates  X,  single  arrivals  in  two  consecutive  slots  contribute  one  third  of  P1(  while 
the  remaining  two  thirds  are  due  to  consecutive  departures  at  the  end  of  a  collision  resolution 
interval.  (2)  As  the  rate  X  of  the  Poisson  input  traffic  increases,  the  interdeparture  distribution 
induced  by  the  algorithm  deviates  further  from  the  Bernoulli  distribution.  In  fact,  as  X  increases, 
the  mass  of  the  interdeparture  distribution  accumulates  at  relatively  small  s  values.  For  example, 

for  X  =  0.4,  we  have  P!  =  0.471  and  £  P,  -  1. 


Remarks  Our  results  showed  that  it  is  generally  wrong  to  conjecture  exponential  interdeparture 
distribution,  (whose  discript  form  is  Bernoulli).  In  fact,  this  distribution  is  far  from  exponential. 
Even  for  small  input  Poisson  rates  X,  the  probability  Pj  does  not  match  the  exponential  fitting. 
We  point  out  that  our  general  approach  in  this  section  and  the  corresponding  regenerative  points 
apply  to  other  algorithms  as  well,  including  the  Capetanakis  dynamic  algorithm.  However,  the 
development  of  the  appropriate  recursions  is  then  an  exceedingly  complex  task.  The  simple 
operations  of  the  proposed  algorithm  present  a  remarkable  advantage,  which  does  not  character¬ 
ize  other  existing  algorithms. 


VI.  Conclusions 

We  presented  a  simple  window  random  access  algorithm  for  systems  with  binary,  collision 
versus  noncollision,  feedback.  We  analyzed  the  algorithm  in  the  presence  of  the  limit  Poisson 
user  model,  and  for  both  its  full  sensing  and  limited  sensing  implementations.  In  addition  to  the 
throughput  and  the  delay  analyses,  we  studied  the  effect  of  feedback  errors  on  the  throughput  of 
the  algorithm  and  the  output  traffic  nterdeparture  distribution,  both  in  the  full  sensing  environ¬ 
ment.  As  compared  to  the  Capetanakis  dynamic  algorithm,  the  proposed  algorithm  is  superior  in 
terms  of  delays  and  insensitivity  to  feedback  errors.  In  contrast  to  the  former,  the  algorithm  can 
also  be  easily  adapted  for  implementation  in  the  limited  sensing  environment,  it  allows  for 
analytical  studying  of  the  output  traffic  interdeparture  distribution,  and  can  be  easily  analyzed 
when  strict  delay  limitations  are  imposed,  [3]. 
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Appendix 

We  first  provide  some  definitions. 


lkjn- 


nk,s: 

fid* 

^d.s* 


P(k,/,5ld): 


Pk(/): 


Given  k  packets  with  counter  values  equal  to  1  and  m  packets  with  counter  values 
equal  to  2,  the  number  of  slots  needed  by  the  algorithm  until  the  first  successful 
transmission,  (and  including  it),  after  the  k-multiplicity  collision  has  been 
observed. 

Given  a  collision  resolution  interval  which  starts  with  a  k-multiplicity  collision, 
the  number  of  length  s  interdeparture  intervals  within  it.  The  length  from  the  ini¬ 
tial  collision  to  the  first  successful  transmission  is  included  in  the  counting. 

Starting  with  a  CRP  at  which  the  lag  equals  d,  d>l,  and  which  follows  a  success¬ 
ful  transmission,  the  number  of  slots  needed  by  the  algorithm  to  reach  the  first 
lag-one  CRP  which  follows  a  slot  containing  a  successful  transmission. 

Starting  with  a  CRP  at  which  the  lag  equals  d,  d>l,  and  which  follows  a  success¬ 
ful  transmission,  the  number  of  length  s  interdeparture  intervals  until  the  first 
lag-one  CRP  which  follows  a  slot  containing  a  successful  transmission.  The  dis¬ 
tance  from  the  initial  CRP  to  the  first  successful  transmission  is  included  in  the 
counting. 

Given  an  arrival  interval  of  length  d,  the  probability  that  there  are  k  arrivals  in  it, 
that  /k  0=8,  and  that  it  takes  l  slots  for  its  resolution,  including  the  initial  collision 
slot. 

Given  a  k-multiplicity  initial  collision,  the  probability  that  it  takes 
/  slots  for  its  resolution,  including  the  initial  collision  slot. 


The  above  definitions  are  needed  for  the  derivation  of  recursions  that  are  pertinent  to  the 
infinite-dimensionality  systems  associated  with  the  quantities  in  (14).  We  first  note  that: 

H  =  E{hj)  ,  C  =  XH  (A.l) 


n=l 


(A.2) 


Auxiliary  Recursions 


The  operations  of  the  algorithm  induce  the  following  recursions: 


(I) 


/l  in  =  0;Vm  ,  P(/k,m=0)=0;  Vk>2,  Vm 


J  1  ,  if  m=l 

^O.m  —  l^m.O  >  P(/0.m-1)-l(),ifm*l 


© 


^kjn  —  1  ^  wn+k- 


;_j  ;  with  probability  j^j  2  k,  k>2 


1  ,  if  k=l  and  s=0 

?0kjn  =  S)  ='  P(/m,o  =  S-l)  ,  if  k=0,  S>1 

r*i  [*]  P  m*k-l  =  s-l).  if  k22,  sal 

i— fi  >  ' 


Jl  ,ifs=l 

ni-s=  io.ifs^l 


f  nk_!  s  ,  with  probability  P(/^0?ts-l) 
*^2  ,  nk  s-  -j  j  +  ^  ^  t  with  probability  P(/^0  =  s-l) 


Nk.s2E{nM}= 


£  P(/i0  =  s-l)  ,  if  s>l 

i=2 

l  +  lP(/i,o  =  OH  if  s=l 

i=2 


'l  , if  /  =  1 

k=0,l ;  Pk(/)=  Q  f  otherwise  ,  P2(/)  =  P(/2.0  =  1-2),  for  />3 


k£2  1  /-  k-i 

/>k+2  ;pk(0=  I  P(/k,0  =  s)Pk_1(/-s-l) 


Given  Poisson  intensity  X, 


P(k,/,p  Id)  =  e_Xd  P(/k.0  =  p)Pk_!  (/-p-1) 


nM  +  m/4*w-P  e  “in- ^(0  ;  te2 


d£A;  m^5  = 


nk-i.* +  m/4  ;  w.p.  P(0 1  d)  n£0  P"(0 1  l)P(k,/  ,p  1 1) ;  k£2 

n+p+2*s 

yk  Z 

1  +  nk-i,s  +  m/<s ;  w.p.  e_X(d+1)—  n£0  eT^PU ^0=P)pk-[0  — p— l);k>2 

*•  n+p+2=s 

.  1 ;  w.p.  Xe~X(d+s_1)  ;  s£2 


For  d  >  A: 


md.*  =  nt5 + md-A+/,s ;  w  P-e  u  —4^-  Pk(/)  ;  tei 

k! 


“  nk-l,$  +  md-n(A-l>f/4  *  W'P-  Z 

p#s-n-l 


e-XA(n+l)  Ml  P(/k  ^p)  Pk_l(/-p-l),  k>l,  if  L  |^-J  > 


=  1  +  nk-i.s  +  md-n(A-iw.s  *  w.p.  Z  e~XA<n+1>-^-  P(f t0  =  s-n-1 )  Pk_!  (/- s+n) ,  k>  i . 


"k-1.5  +  mi.» ;  w.p.  e  ± - - - J- 


X  P(/k,o  =  P)Pk.i(f-p-l);k>2 


=  1  +nk-l4  +  tn/4 


-4h^1  d_A] 

;  w.p.  e  - - - - —  P  /k.o  =  s-l-L^J 

»  - 

Pk-.  :k22 

k  « 


nk-l4  +  m/,»  *  w-P-e 


-^rp+1) 

A_1  n£0 


e->«^.p(/t0  =  p)Pk_,(/-p-1);k22 


n+2+p4-^^j  « 


—A,  d-»{  — — — j  +1 

=  H-n^  +  m^  ;w.p.e  A_1  I  e'^  P  /w  =  s-n-2-j.  . 

n>n  A— I 


.Pk_!  /+n-s4{^-^+l  ;k>2 
A-l 


=  1  ;  w.p.Xc_x<d+*_1) ;  if  s-2-J.  -2-^j  £0 

A-l 


JA  1  -*■  J  A 

=  1  ;  w.p.X  |  (A-l)  e  ;  if  s  =  l-t{ 

A— 1  A— 1 


Let  us  define. 


A  i x£0 

u<«=  «,,<0 

p5<'>=  I  «-*  W 

kit  K‘ 


N8.s=  I  Nk,s  e~*  ,  Ng  i  =  1-e-® 

k2l  K! 


ps,=  I e-t  |r  p<'k.o= 8) 


(A.2) 


Then,  using  the  above  defined  quantities,  and  the  recursions  on  md  s,  we  easily  find: 


For  d  £  A: 


4  —  AX1 

NU*+  I  Px  (/)!-*- 


>*  v  VV 


'  vWV 


p 

nl: 


rti 


+  — T  [NU  "  *\s-ij  +  U(s-2)e--|e^*-^  £  e™  P„m  + 


+  Xc~Ms~l)  -  \-±- 


(A.3) 


For  d>A  and  p  ^  [ 'T-Tj  ,  p  =  0,l, ...: 

-  A— 1 


^d,s  ^aAs  2  M<1 -A+/^Pxa(^) 
121 


+  e-XAp-|  N 


X(d+p-pA),s  f*X<d+p-pA),s-l  X  M,,  Px(d+p-pA)(0 

123 


_-X(d+p)  f 

J-  “ — —]  -  P^.!  +  £  M,tS  PX  (/) 

1-C  1  123 


UCp-J  e^-(1^eJAP)  [n^-P^J  +  £  e-Un£Md_n(A_1)+/iSPXA(/)l 

L  1-e  L  J  lSSnSp  12 1  J 


+  UCp-DUCs^c"^'1)  £  e^P^, 

m=s-l-min(p,  s-1) 


-  (d+p-pA)eXpA  + 


+  UCs-l-pJc-^P^^!^ 


+  U(s-2-p)  c_X(d+s“2)  *  £  P  P^ 

m=0 


(A.4) 


Bounds 


For  the  numbers  N^,,  we  used  the  following  bounds: 


0  S  Nk  s  <,  k-1  ;  Vs 

Regarding  the  numbers  Hd,  we  used  the  methodology  in  [4],  and  proved  that, 


W.VTji’ji  V. 
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cqd+p,  ^Hd^ctud  +  p,,  ,  d£l 


(A.6) 


where. 


cli  =au  =  [A-E{/ 1  A}]-1  E{/ 1  A} 


l<d<A 


Q(d)  ,  pu  =  max[^-au,  1J^AQ(d)j 


Q(d)=  [Xde-Xdj  \  E(/  Id}  +a^E{/  Id}  -  d-Xde_Xxl 


Bounds  on  the  numbers  Md  s  can  be  developed  similarly  with  those  for  the  numbers  Hj.  The 
former  are  significantly  more  complicated,  however.  Instead,  we  used  the  following  simpler  and 
intuitively  clear  bounds,  where  Hd  denotes  the  upper  bound  on  the  quantity 


0  ^  M^s  <Hd 


(A.7) 


We  used  the  bounds  in  (A.7),  for  d  £  30. 
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Window  Selection  in  the  Limited  Sensing  Environment 
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Figure  3 


Expected  Delays  Induced  by  the  Proposed  Algorithm 
in  the  Full  Feedback  Sensing  and  the  Limited  Feedback 
Sensing  Environments 
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